ANSWERS (OPEN-CHANNEL FLOW NOTES) AUTUMN 2023 


Section 1.2 


Example. 

The discharge in a channel with bottom width 3 m is 12 m° s!. If Manning’s n is 0.013 m7 s 
and the streamwise slope is 1 in 200, find the normal depth if: 

(a) the channel has vertical sides (i.e. rectangular channel); 

(b) the channel is trapezoidal with side slopes 2H:1V. 


b = 3 m (base width) 


Q = 12 m? 54 
n = 0.013 m71/3 s 
S = 0.005 
(a) 
Discharge: 
Q=VA 
where, in normal flow, 
Vio bh h 
= — /2 A= bh, R, = —— = — 
ra NO t= 2h ~ 1+2h/b 
Hence, 
5/3 
1 bh” 51/2 


~ n(1 + 2h/b)2/3 
Rearranging as an iterative formula for h: 


Here, with lengths in metres, 
h = 0.8316 (1 + 2h/3)?/5 
Iteration (from, e.g., h = 0.8316) gives 


hy = 1.024 m 


Answer: 1.02 m. 


(b) Geometry: trapezoidal cross-section with base width b, surface width b + 2 x (2h) and two 


sloping side lengths 4h? + (2h)? = hv5. 
Area and wetted perimeter: 
A= =(b +b+4h)h =h(b+2h) =hb(1+2h/b) 


P=b+2hv5 
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Hydraulic radius: 


nQ 


ae A 
np 
Discharge: 
Q=VA 
Hence, 
1 
Q=- 
n 
nQ 
=> — Z 
bVS 
=. Qa (= 
bVS 


120 1+ 2h/b 


5/3 


o h(b+2h) ( 1+ 2h/b 
b+2hv5 1+ 2V5h/b 
== RI S24 
n 


2/3 
S1/2hb(1 + 2h/b 
1+ a) ( /b) 


(1 + 2h/b)5/? 
(1 + 2V5h/b)2/3 


i (1 + 2V5h/b)?/° 


1+2h/b 


Here, with lengths in metres, 


h = 0.8316 


(1 + 1.491h)?/° 
1+ 2h/3 


Iteration (from, e.g., h = 0.8316) gives 
h, = 0.7487 m 


Answer: 0.749 m. 
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Section 1.4 


Example. 

The discharge in a rectangular channel of width 6 m with Manning’s n = 0.012 m 
24 m° s™!. If the streamwise slope is 1 in 200 find: 

(a) the normal depth; 

(b) the Froude number at the normal depth; 

(c) the critical depth. 


“1/3 ¢ is 


State whether the normal flow is subcritical or supercritical. 


b=6m 

n = 0.012 m71 s 
Q = 24m? s7! 

S = 0.005 


(a) 
Discharge: 
Q=VA 
where, in normal flow, 
Lati bh h 
= [2 A = bh, R = —_ = —— 
i o h b+2h  1+2h/b 


Hence, 
1 bhè”? 
SEE o 
n (1 + 2h/b)2/3 


or, rearranging as an iterative formula for h: 


_ ("Q R 2/5 
h= a. (1 + 2h/b) 


Here, with lengths in metres, 


h = 0.7926 (1 + h/3)2/5 


1/2 


Iteration (from, e.g., h = 0.7926) gives 
h, = 0.8783 m 


Answer: 0.878 m. 


(b) At the normal depth, h = 0.8783 m: 
Q _ 24 
A 6X 0.8783 


V 4,554 
Fr = == Ss ————._ = 1551 
Vgh v9.81 x 0.8783 


V= = 4,554 m s71 


Answer: 1.55. 
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(c) The critical depth is that depth (at the given flow rate) for which Fr = 1. It is not normal 
flow, and does not depend on the slope S or the roughness n. 


V 
Fr = —— 
Vgh 
where 
Q 
V=% in general 
7 Gng ) 
Q q S ENE 
or V = oh R (for a rectangular channel; q is the flow per unit width) 


Hence, for a rectangular channel, 


Fr2 = (q/h)’ = a 
gh gh? 


For critical flow, Fr = 1 and so 


2\ 1/3 
E 
g 


Here, the flow per unit width is q = 24/6 = 4 m? s“1, so that 


42 \1/3 
he = (a) = 1.177 m 


Answer: 1.18 m. 


The normal depth is supercritical because, when h = h,,, then Fr > 1 (part (b)). 


Alternatively (and often more conveniently), the normal depth here is supercritical because 
hy < he; so speed V is larger, and depth h is smaller in normal flow than critical flow, so that 


Fr =V/./gh must be greater than 1. 
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Section 2.2 


Example. 

A 3-m wide channel carries a total discharge of 12 m 
(a) the critical depth; 

(b) the minimum specific energy; 

(c) the alternate depths when E = 4 m. 


3 s}, Calculate: 


b=3m 
Q =12m? s7! 
(a) 
Discharge per unit width: 
12 
q= 7 = z =4 m2 s71 


Then, for a rectangular channel: 


ee Tee. 24497 
=y ~ (981 a 


Answer: 1.18 m. 


(b) For a rectangular channel, 


3 
he =3X1177 =1.766m 


Answer: 1.77 m. 


(c) As E > E¢, there are two possible depths for a given specific energy. 


Ve Q q 
E=h+— where V=% =- (for a rectangular channel) 
2g A h 
2 
q 
E=h 
=> + TI 
Substituting values in metre-second units: 
0.8155 


For the subcritical (slow, deep) solution, the first term, associated with potential energy, 
dominates, so rearrange as: 
0.8155 


RRAS 
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Iteration (from, e.g., h =4) gives h = 3.948 m. 


For the supercritical (fast, shallow) solution, the second term, associated with kinetic energy, 
dominates, so rearrange as: 


a 0.8155 
~ | 4—h 


Iteration (from, e.g., h = 0) gives h = 0.4814 m. 


Answer: 3.95 m and 0.481 m. 
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Section 2.3.1 


Example. (Exam 2020) 
(a) Define: 
(i) specific energy 
(11) Froude number 
for open-channel flow. What is special about these quantities in critical conditions? 


A long, wide channel has a slope of 1:1000, a Manning’s n of 0.015 m9 s and a discharge of 

3 m° s™! per metre width. 

(b) Calculate the normal and critical depths. 

(c) In a region of the channel the bed is raised by a height of 0.8 m over a length sufficient 
for the flow to be parallel to the bed over this length. Determine the depths upstream, 
downstream and over the raised bed, ignoring frictional losses. Sketch the key features 
of the flow, indicating all hydraulic transitions caused by the bed rise. 

(d) In the same channel, the bed is lowered by 0.8 m from its original level. Determine the 
depths upstream, downstream and over the lowered bed, ignoring frictional losses. 


Sketch the flow. 
(a) 
(i) Specific energy is the head relative to the local bed of the channel: 
V2 
EF=h+=— 
+ Zg 


Gii) The Froude number is 


Fr = —— 


In critical conditions Fr = 1 and the specific energy is the minimum for that discharge. 


(b) 

S = 0.001 

n = 0.015m?/35 
q = 3m? s* 


Normal depth: 
Discharge per unit width: 


1 
q=Vh, where V =-= Re 351/2 (Manning), Rp =h (“wide” channel) 
n 
1 
> q=—h/351/2h 
n 


he/3./5 
7 n 


> og 
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3/5 


L a -i ens x *) T 
= f — = f; m 
” VS V0.001 


Critical depth: 


he = AS (2y = 0.9717 
ar ~ \o81 T Š 


Answer: normal depth = 1.24 m; critical depth = 0.972 m. 


(c) To determine the type of behaviour over the raised bed, compare the total head under critical 
conditions (the minimum energy necessary to get over the weir at this flow rate) with that 
available in the approach flow. 

Critical 


he = 0.9717 m 
3 
E,=5he =1458 m 


Zp = 0.8 m 


He =Z +E, = 2.258 m 


Approach Flow 
Because the channel is described as “long” it will have sufficient fetch to develop normal flow; 
hence the approach-flow head is that for the normal depth (h = 1.236 m): 


V2 
2g 


Ha = Ea = hn t 
q? 
= h, += 
” 2ghê 
32 
= 1.236 + —_____ 
t 2 x 9.81 x 1.2362 
= 1.536 m 


At the normal depth the available head (H4) is less than the minimum required to get over the 
bed rise (H¿). Hence the water depth must increase upstream (“back up”), to raise the head 
immediately upstream. Thus: 

° critical conditions do occur; 

e the total head in the vicinity is the critical head (H = H, = 2.258 m). 


Over the raised bed there is a hydraulic transition, so the depth over this is critical: h = hę = 
0.9717 m. 


Just up- or downstream, 
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2 


Pfam h ee 
=F= — where =— 
2g h 
2 
q 
H=h 
=> tagh 
0.4587 


Upstream, rearrange for the deep, subcritical, solution: 
0.4587 

h2 
Iteration (from, e.g., h = 2.258) gives h = 2.160 m. 


h = 2.258 — 


Downstream, rearrange for the shallow, supercritical solution: 


T 0.4587 
~ 12.258- h 


Iteration (from, e.g., h = 0) gives h = 0.5127 m. 


normal GVF RVF 


normal 


Since the preferred (i.e. normal) depth is subcritical, there must be a downstream hydraulic 
jump. (A quick calculation shows that the upstream depth for this jump is greater than hz, so 
there is indeed a length of GVF between the area of bed rise and the jump.) 


Answer: depths upstream, over, downstream of the raised bed: 2.16 m, 0.972 m, 0.513 m. 


(d) The flow does not require additional energy to pass a depressed section; hence, the total 
head throughout is that supplied by the approach flow (H = Hg = 1.536 m) and the flow 
remains subcritical. The depths just upstream and downstream of the lowered section are those 
in the approach flow; i.e. normal depth. 


As bed height z, decreases, specific energy E must increase to maintain the same total head. 
In the lowered section: 


H = Zb +E 
=> 1.536 = —0.8 + E 


=> E = 2.336 m 
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Then 
2 


Pein h v=) 
= — where =— 
2g h 
2 
_ q 
> DERT ae 


0.4587 
As we require the subcritical solution, rearrange as 
0.4587 
h2 
Iteration (from, e.g., h = 2.336) gives h = 2.245 m. 


h = 2.336 — 


(Note that this is the depth of the water column. The actual surface level here is 
Zs = —0.8 +h = 1.445 m 


so the overall water level also rises in this section.) 


normal RVF 


normal 


Answer: depths upstream, within, downstream of the lowered section: 1.24 m, 2.24 m, 1.24 m. 
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Section 2.3.1 


Example. 

A long channel of rectangular cross-section with width 3.5 m and streamwise slope 1 in 800 
carries a discharge of 15 m° s™!. Manning’s n may be taken as 0.016 m~™!® s. A broad-crested 
weir of height 0.7 m is constructed at the centre of the channel. Determine: 

(a) the depth far upstream of the weir; 

(b) the depth just upstream of the weir; 

(c) whether or not a region of supercritical gradually-varied flow exists downstream of the 


weir. 
b=3.5m 

S = 0.00125 
Q=15:m*-s7* 


n=0.016m7!/3 s 
Zweir = 0.7 m 


(a) The depth far upstream is normal since the channel is described as “long”. For normal flow 
in a rectangular channel: 


Q=VA 
where: 
1 2/301 bh h 
ee ores a ae rT 1+ 2h/b 
Hence, 
1 bk” 


1/2 


= ——— zS 
n (1 + 2h/b)2/3 
or, rearranging as an iterative formula for h: 


3/5 


Ke (=) (1 + 2h/b)2/5 


Here, with lengths in metres, 
h = 1.488 (1 + 0.5714h)?/5 
Iteration (from, e.g., h = 1.488) gives 
h = 2.023 m 


Answer: 2.02 m. 


(b) To establish depths near the weir we need to know the flow behaviour at the weir. Compare 
the energy in the approach flow with that under critical conditions. 


Approach flow 


Hydraulics 3 Answers (Open-Channel Flow Notes) - 11 Dr David Apsley 


h = 2.023m (from part (a)) 
Q Q 15 


“a bh a 
Specific energy in the approach flow: 
v? 2.118? 
ar = 2.023 + 981 = 2.252 m 


Referring heads to the undisturbed bed near the weir: 
Ha = Ea = 2.252 m 


Critical conditions 


2\ 1/3 15 
--(5) q= =77 =4286m? s- 


E.==h, =Żł2x1.233 =1.850m 


NIW 


He = Zweir t Ee =0.7+1.850 =2.550m 


Since the head required to flow over the weir (H; = 2.550 m) exceeds that in the approach 
flow (Hg = 2.252 m), the depth just upstream of the weir must increase and the flow back up. 
The total head at any position in the vicinity of the weir is H = H; = 2.550 m. 


Just upstream and downstream of the weir (i.e. at undisturbed bed level): 


H=E rea V= Q 
2g ~ bh 
a Q? 
= EP pme 
0.9362 


The depth just upstream is the deep, subcritical solution. Hence, rearrange as 


0.9362 
h = 2.550 = “hR 
Iteration (from, e.g., h = 2.550) gives 
h = 2.385 m 


Answer: 2.39 m. 
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(c) Since the normal flow is subcritical, the flow must return to it via a hydraulic jump on the 
downstream side of the weir. 


If the flow in the vicinity of the weir is unaffected by the hydraulic jump the flow goes smoothly 
supercritical on the downstream side, with total head H = 2.550 m (equation (*)). Rearranging 
to get an iterative formula for the supercritical solution: 


_ | 0.9362 
~ 12.550—h 


Iteration (from, e.g., h = 0) gives 
h = 0.7141 m 


Denote by subscripts A and B respectively the conditions upstream and downstream of the 
hydraulic jump. On the downstream side conditions may be assumed normal, since the channel 
is “long” and hence there is sufficient fetch to develop the preferred depth: 


hp = 2.023 m 
Ve = 2.118 ms" (from part (b)) 
Vs 2.118 
foes Sa S751 


«ghg v9.81 x 2.023 


Hence, from the hydraulic-jump relation for the sequent depths: 


h 2.0 
ha == (-1+ 1 + 8Fr2) =e 1+/1+8x 0.47542) = 0.6835 m 


Any gradually-varied supercritical flow downstream of the weir would increase in depth until 
a hydraulic jump occurred (see the lectures on GVF). Since the depth downstream of the weir 
is already greater than any sequent depth upstream of the hydraulic jump, no such increasing- 
depth GVF is possible and the hydraulic jump must actually occur at (or just before) the 
downstream end of the weir. 
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Section 2.3.1 


Example. 

A reservoir discharge is controlled by a weir of width 8 m and discharge coefficient 0.9. 

(a) Calculate the flow rate over the weir when the freeboard is 0.65 m. 

(b) Assuming negligible inflow and a constant plan area of 1.5 km? for the reservoir, 
calculate the time in hours taken to reduce the level of the reservoir by 0.4 m. 


(a) The weir controls the flow, so there must be critical flow over the weir and the freeboard 
corresponds to critical head: 


he = ho 


2 
> q= Ja(5ro) 


For the full width b = 8 m, and with a discharge coefficient cg = 0.9, 


2 3 
O=csqb = tab |g (Éh) = 12.28n2/? 


When hy = 0.65 m, 
Q = 6.435 m? s71 


Answer: 6.44 m? s™!. 


(b) By continuity: 
dV 
dt = Qin — Qout 


Here, dV = Aws dh, Qin = 0 and Qout can be written as a function of h (dropping subscript 0) 
as above. Hence, 


h 
Aws zz = 712-28 h3/2 
Writing Aws = 1.5 km? = 1.5 x 10° m? and separating variables: 
1.5 x 10° 
—_ 3/2 — 
12.28 hœ” -dh = dt 
Integrating between t = 0 (where h = 0.65 m) and t = T (where h = 0.25 m): 
0.25 T 
-122100 Í h-3/2dh = | dt 
0.65 0 
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=D 
244200 (—— pee ) T 
> — ramal 
v0.25 0.65 


=> T = 185500s = 51.53 hr 


Answer: 51.5 hr. 
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Section 2.3.2 


Example. 


A venturi flume is placed near the middle of a long rectangular channel with Manning’s n = 
0.012 m~1/3 s. The channel has a width of 5 m, a discharge of 12.5 m? s! and a slope of 
1:2500. 


(a) Determine the critical depth and the normal depth in the main channel. 

(b) Determine the venturi flume width which will just make the flow critical at the 
contraction. 

(c) If the contraction width is 2 m find the depths just upstream, downstream and at the 
throat of the venturi flume (neglecting friction in this short section). 

(d) Sketch the surface profile. 


n=0.012m-'/3 s 
b = 5 m (main channel) 
Q = 12.5 m? s71 


S$=4x10%* 
(a) 
In the main channel, 
12.5 
q= irs = 2.5 m? s7! 


Normal Depth 


Q=VA 
where: 
1 2/304 bh h 
aa ae h b+2h  1+2h/b 
Hence, 
1 bh?” 


1/2 


= —— sS 
n (1 + 2h/b)2/3 
or, rearranging as an iterative formula for h: 


_ ("Q G 2/5 
h= =a (1 + 2h/b) 


Here, with lengths in metres, 
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h = 1.275 (1 + 0.4h)?/5 
Iteration (from, e.g., h = 1.275) gives 
hp = 1.546 m 


Answer: critical depth = 0.860 m; normal depth = 1.55 m. 


(b) The flow will just go critical if the head in the throat (H,) is exactly equal to that in the 
approach flow (H,). Measure heads relative to the bed of the channel in the vicinity of the 
venturi. 


Critical Head 
/3 
3 3 qV 0 
H, = E, = she -3(2) r where dm = — 


(Note that the critical depth is different at the throat to that in the main channel, due to the 
narrower width.) 


P 3 Q? 1/3 
~ ZNgbR, 


Approach Flow 
The approach flow is normal, since the channel is “long”. Hence, 
ha = 1.546 m 
Q 12.5 
VW =—— == =1.617ms™! 
ac ph - 5X1546 tee 
Ha =E, =h Va = 1.546 + soe = 1.679 
BEE. TES ag. A IROL eee 


For the flow just to go critical at the throat, 


He = Ha 
3 / 02 \¥3 
=>. (3) = 1.679 
2\gbin 
12.5 
Š bm = ry, = 3.370 m 


é x 1.679) V381 


Answer: 3.37 m. 
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(c) If the throat width is reduced further, then the flow will back up and undergo a critical 
transition at the throat. 


At the throat, 


12.5 
fate =—— =6.25m’s"! 
bm 
3 (q2,\"" 6.252\ 
he=3(@) - (Sr) = 1.585 m 


3 3 
Ez = zhe = 7 x 1.585 = 2.378 m 


He =Z, +E, =0+2.378 =2.378m 


The head throughout the venturi will be the critical head (H = H, = 2.378 m). 


Anywhere in the flume, 


y? Q 
H aT where Zs V bh 
= Q? 
> KERETE 


At the throat the depth will be the critical depth there; i.e. h = h, = 1.585 m. 


Just upstream and downstream, b = 5 m; hence, 


0.3186 
2.378 = h + “R 


Upstream 
Rearrange for the deep, subcritical solution: 
0.3186 
h = 2.378 — “R 


Iteration (from, e.g., h = 2.378) gives h = 2.319 m. 


Downstream 
Rearrange for the shallow, supercritical solution: 


he 0.3186 
(2.378 —h 
Iteration (from, e.g., h = 0) gives h = 0.4015 m. 


Answer: depths upstream, in the throat, downstream = 2.32 m, 1.59 m, 0.401 m. 
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WATER PROFILE 


critical 
t 
—> 
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Section 2.3.3 


Example. 

The water depth upstream of a sluice gate is 0.8 m and the depth just downstream (at the vena 
contracta) is 0.2 m. Calculate: 

(a) the discharge per unit width; 

(b) the Froude numbers upstream and downstream. 


hy = 0.8 m 
h> = 0.2 m 


(a) Assuming total head the same on either side of the gate: 


V2 V2 
Zsı + 2g = Zş2 + 2g 
Substituting z, = h and V = q/h: 
q? q? 
h +—> = h, t+ 
1 2gh? °? 2gh? 


From the given data, in metre-second units: 


0.8 + 0.0796q? = 0.2 + 1.2742q? 
> 0.6 = 1.1946q? 


= q = 0.7087 m? s71 


Answer: 0.709 m? s™!. 


(b) Use, on each side of the gate, 


to get 
V, = 0.8859 ms"? 
V, = 3.544 m s™t 


and then 
Fr, = 0.3162 
Fr, = 2.530 


Answer: Froude numbers upstream, downstream = 0.316, 2.53. 
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Section 2.3.3 


Example. 
A sluice gate controls the flow in a channel of width 2 m. If the discharge is 0.5 m? s™! and the 
upstream water depth is 1.5 m, calculate the downstream depth and velocity. 


b=2m 
Q = 0.5 m? s71 
h =1.5m 


Use upstream conditions to get total head. Then, assuming no losses, find the supercritical flow 
with the same head. 


Total head (either side): 


y? Q 
H = — h =h d = — 
Zs + Zg where Zs an V A 
i e a RAO 
= OT gbthe | h? 


The upstream depth h; = 1.5 m gives 


3.186 x 1073 ; 
H=1.5+ pee = 1.501 m (dominated by h,) 
Hence, 
3.186 x 1073 
1.501 = h + BF 


Rearrange for the shallow, supercritical solution: 


ie 3.186 x 10-3 
~ | 1.501 —h 


Iteration (from, e.g., h = 0) gives 
h, = 0.04681 m 


Q 0.5 


ee ee ee -1 
+o. Ixo O Ms 


V2 


Answer: downstream depth = 0.0468 m; velocity = 5.34 m s™!. 
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Section 2.4 


Example. (Exam 2018) 


Water flows at 0.8 m° s™! per metre width down a long, wide spillway of slope 1 in 30 onto a 


wide apron of slope 1 in 1000. Manning’s roughness coefficient n = 0.014 m 


slopes. 


-1/3 s on both 


(a) Find the normal depths in both sections and show that normal flow is supercritical on 


the spillway and subcritical on the apron. 


(b) Baffle blocks are placed a short distance downstream of the slope transition to provoke 
a hydraulic jump. Assuming that flow is normal on both the spillway and downstream 
of the hydraulic jump, calculate the force per metre width of channel that the blocks 


must impart. 


(c) Find the head loss across the blocks. 


Sı = 1/30; S, = 1/1000 
q = 0.8 m? s71 
n=0.014m71/3 s 


(a) Normal flow: 


1 ; 
q=Vh, where V =- Re! $51/2 (Manning), R,=h (“wide” channel) 
n 


=> q = 1 2/351/2h 
n 


5/3 
J e h an 
ngq\3/5 
= t=) 
For the two slopes this gives 
h, = 0.1874 m 
h, = 0.5365 m 


Answer: depth on spillway = 0.187 m; depth on apron = 0.536 m. 


For subcritical/supercritical: 
Method 1 (use critical depth) 
The critical depth is 


h. = T N ps m = 0.4026 
e= Ng ~ \ogt eee 
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On the spillway, h, < he: this is shallower than critical flow (where Fr = 1) and hence faster; 
both ensure Fr; > 1, so supercritical. 

On the apron, hz > họ: this is deeper than critical flow and hence slower; both ensure Fr, < 1, 
so subcritical. 


Method 2 (find Froude numbers) 


V 
foe: oi 


Jah ~ Jar 


Applying this for both depths we find Fr, = 3.148 (supercritical) and Fr, = 0.650 
(subcritical). 


(b) 


The corresponding velocities are deduced from V = q/h, whence: 
V, = 4.269 mst 
V, = 1.491 m s™t 


Let f be the magnitude of the force per unit width exerted by the fluid on the blocks and, by 
reaction, the blocks on the fluid, which is clearly acts in the upstream direction. 


On each side of the blocks the hydrostatic pressure force is given by 


average pressure X area 


1 
or J pghxh (per unit width) 


Hence, from the steady-state momentum principle: 


force = rate of change of momentum (mass flux x change in velocity) 
1 2_1 2 
-f +5Ppghi — 5 pgh2 = pa = V1) 
Hence, 
1 2 2 
f =3P9 Chi — hz) + pqs — V2) 


= —1240 + 2222 
= 982 N m“! 


Answer: 982 N m™!. 


(c) Where hydrostatic, head in open-channel flow is given by 
2 


H =z; + — 
Zs 2g 


Here, relative to the bed of the apron, z, = h. Hence, 
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head loss = h, — h Pie ibe 
ead loss = h; — h, 2g 
= —0.3491 + 0.8156 
= 0.4665 m 


Answer: 0.467 m. 
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Section 2.4 


Example. 
A downward step of height 0.5 m causes a hydraulic 
jump in a wide channel when the depth and velocity 
of the flow upstream are 0.5 m and 10 ms", 
respectively. 


(a) Find the downstream depth. 


(b) Find the head lost in the jump. 


(a) The downstream depth can be deduced from the momentum principle if the reaction force 
from the step is known. The approximation is that this is the same as would occur if it were in 
equilibrium with a hydrostatic pressure distribution here. 


Flow rate per unit width: 
q=Vıh =10x0.5 =5 m?s“! 


Steady-state momentum principle 
force = rate of change of momentum (mass flux x change in velocity) 


Per unit width: 


1 2 1 2 
z PI Ch +4) — z Pgh: = pq(V2 — V1) 
Since 


DE 
oh, h 


this gives, in metre-second units, 


V2 


5 
4905(1 — h) = 5000(-— — 10) 
2 


5 
> 1-—hA2 =1.019(—- 10) 
hy 


> 1119=h24+ 


Rearrange for the deep, subcritical solution: 


5.095 
hy = [11.19- 
hy 


Iterating (from, e.g., hy = V11.19) gives hz = 3.089 m. 


Answer: 3.09 m. 
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(b) Head either side is given by: 
2 
29 
The datum is not important as it is only the difference in head that is required. For convenience, 
measure Z relative to the bed of the expanded part. Then, 


H = Z; + 


Zs, = 1m (note: water surface level, not depth), V = 10 ms7™?! = H = 6.097 m 


Zs. = 3.089 m, V> = q/h, = 1.619ms7! => H, = 3.223 m 
Hence, 


head lost = 6.097 — 3.223 = 2.874 m 


Answer: 2.87 m. 
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Section 3.6.2 


Example. 

A long, wide channel has a slope of 1:2747 with a Manning’s n of 0.015 m™® s. It carries a 
discharge of 2.5 m? s™! per metre width, and there is a free overfall at the downstream end. An 
undershot sluice is placed a certain distance upstream of the free overfall which determines the 
nature of the flow between sluice and overfall. The depth just downstream of the sluice is 0.5 m. 


(a) Determine the critical depth and normal depth. 


(b) Sketch, with explanation, the two possible gradually-varied flows between sluice and 
overfall. 


(c) Calculate the particular distance between sluice and overfall which determines the 
boundary between these two flows. Use one step in the gradually-varied-flow equation. 


So = 1/2747 = 3.640 x 1074 
n = 0.015 m7!/3 s 
q = 2.5 m? s71 


(a) 
Critical Depth 


Normal Depth 
1 : 
q=Vh, where V=- Re ss 2 (Manning), Rp „=h (“wide” channel) 
n 


1 
=  q=zh S h 


n 


nq ) (ons x =) T 
E —} ee __ = . m 


oS (a 1/2747 


Answer: critical depth = 0.860 m; normal depth = 1.50 m. 


(b) The depth just downstream of the sluice is supercritical (0.5m < h,). However, the 
preferred depth is subcritical (h, > he). Hence, if the channel is long enough then there will 
be a downstream hydraulic jump, with the flow depth then decreasing to pass through critical 
again near the overfall. 
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If the channel is too short, however, the region of supercritical flow from the sluice will extend 
to the overfall. 


hydraulic jump 


normal 


supercritical — — critical 


— — — — normal 


supercritical — — — — critical 


TN 


(c) As the channel shortens, the depth change across the hydraulic jump diminishes. The 
boundary between the two possible flow behaviours occurs when the supercritical GVF just 
reaches critical depth at the overfall (i.e. the limiting depth change across the hydraulic jump 
is zero). 


As the flow is supercritical, integrate the GVF equation forward from the downstream side of 
the sluice gate (where h = 0.5 m) to the overfall (where h = h, = 0.8605 m). Use 1 step. 


GVF equation: 
dh = So — Ss 
dx 1—Fr? 
For the direct-step method invert the GVF equation: 
dx 1-— Fr? dx 
T T and Ax = (=) Ah 


For the working, write the derivative as a function of h; (all lengths in metres). 


V q q? 0.6371 


Fr = —— = — 2ean E aoe 
Jgh Jgh? ae h3 
nq\2 1.406 x 1073 
s= (=z) = h10/3 


Ah = 0.8605 — 0.5 = 0.3605 


Hydraulics 3 Answers (Open-Channel Flow Notes) - 28 Dr David Apsley 


Working formulae: 
A (=) Ah 
x = |— 
dh mid 
where 
0.6371 
ee | Ah = 0.3605 
dh -, 1406x107? ied 
, dx 
l hi Xi hmia ass Ax 
0 0.5 0 
0.6803 217.1 78.26 
1 0.8605 78.26 


Answer: 78 m. 
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Example. (Exam 2022) 

A long rectangular channel of width 2.5 m, slope 0.004 and Manning’s roughness coefficient 
n = 0.022 m~t s carries water at 4m? s71. Temporary works narrow the channel at one 
location to 1.1 m for a short distance. 


(a) Find the normal depth in the main channel and show that the slope is hydraulically mild. 


(b) Show that a hydraulic transition takes place at the narrow point and find the depth just 
downstream of the narrowed section, confirming that supercritical flow is possible here. 


(c) Use two steps in the gradually-varied-flow equation to estimate the distance from the 
end of the narrow section to the downstream hydraulic jump. 


(a) 
b=2.5m (bm = 1.1m in the narrow section) 
So = 0.004 
n = 0.022 m71 s 
Q = 4m? s™t 
Discharge: 
Q=VA 
where, in normal flow: 
ee ee bh h 
z= — /2 = = —— = — 
ee et AEO Hae 1+ 2h/b 
Hence, 
i, “phere 


== SL /2 
n (1 + 2h/b)2/3 


or, rearranging as an iterative formula for h: 


3/5 


= ng 2/5 * 
h= (5) (1 + 2h/b) (*) 


When S = Sọ = 0.004, substitution of numerical values in (*) yields an iterative formula 
h = 0.7036(1 + 0.8h)?/° 

Iteration (from, e.g., h = 0.7036) gives 
h, = 0.8690 m 


At the normal depth the velocity is V, = Q/(bh,) = 1.841 ms? and the Froude number is 
Vn 
Ihn 


This is less than 1, i.e. subcritical. Hence, the slope is hydraulically mild. 


= 0.6305 


Fr = 


Answer: 0.869 m. 
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(b) At the throat of the narrow section the critical depth is 


1/3 
ma (=) 
€ g 


where the flow rate per unit width qm = Q/bm = 3.636 m? s~+. Hence 
he = 1.105 m 


The critical head is 


3 
He =%y +E: =O+5h, =1.658m 


The approach flow head is based on normal flow (since the channel is long): 


2 
n 


V, 
Ha = In +37 = 1.042m 


This is less than the critical head (the minimum energy required to pass this flow through the 
constricted section). Hence, upstream of the narrow section the flow must back up, increasing 
in depth to provide the required energy, and so a hydraulic transition takes place. The total head 
in the vicinity of the narrowed section is 


H=H, =1.658m 


Downstream of the constricted section we require supercritical flow at the same head: 
Q? 
+ ————— 
2gb?h? 


i.e., in m-s units, 


H=h 


0.1305 
1.658 =h+—, 


Rearranging for the supercritical (shallow) solution: 


ps 0.1305 
~ 11.658—h 


Iteration (from, e.g., h = 0) gives 
h, = 0.3113 m 


To show that this supercritical flow actually occurs we must show that it is less than the depth 
upstream of the hydraulic jump (since supercritical GVF on a mild slope increases in depth 
with distance). As the channel is “long”, the hydraulic jump must cause a return to the normal 
depth and the depth on the upstream side of the jump is the sequent depth. With Fr, = 0.6305 
as calculated in part (a), 


h 
hy=—(-1+ 1+8Frå) = 0.4539 m 
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Since hz < h; a depth-increasing length of supercritical flow exists. 


Answer: 0.311 m. 


(c) GVF from just downstream of the narrow section (hstart = 0.3113 m) to the upstream side 
of the hydraulic jump (hena = 0.4539 m). The flow is supercritical so work in the direction of 


the flow. Using two steps the depth increment is 
Rend — Astart _ 0.4539 — 0.3113 


Ah = = 0.0713 m 
Nsteps 2 
GVF equation: 
dh So —S¢ 
dx 1-—Fr? 
For the direct-step method invert the GVF equation: 
dx 1-—Fr? i dx 
dh So- Sp ü ae ea oH 
For the working, write the derivative as a function of h. 
V 2 
EE ES i T _ 0.2610 
q gh bh;/ gh b2? gh? h? 
nQ y? _, (1 + 0.8h)43 
Working formulae: 
A (=) Ah 
x =|— 
dh mid 
where 
0.2610 
eee at Ah = 0.0713 
Tonan nanona Y =0. 
dh 14 — 1.239 CH] x 10-3 
i hi Xi Amia (dx/dh) mia Ax 
0 0.3113 0 
0.34695 96.27 6.864 
1 0.3826 6.864 
0.41825 87.70 6.253 
2 0.4539 13.12 


Answer: 13.1 m. 
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Example. (Exam 2021) 
A long rectangular channel of width 2.2 m, streamwise slope 1:100 and Chézy coefficient 
80 m!” s~! carries a discharge of 4.5 m° s™!. 


(a) Find the normal depth and critical depth and show that the slope is steep at this 
discharge. 


(b) An undershot sluice gate causes a hydraulic transition in this flow. The depth of parallel 
flow downstream of the gate is 0.35 m. Find the depth immediately upstream of the 
gate and sketch the flow. 


(c) Using 2 steps in the gradually-varied-flow equation, find the distance between the gate 


and the hydraulic jump. 
(a) 
b = 2.2 m 
So = 0.01 
C = 80m? s7! 
Q = 4.5 m? s7! 


For the normal depth, 


Q=VA 
where 
bh h 
= 1/2 61/2 z - = 
V = CRE So, BS LES on TFN 
Hence, 
h 1/2 ‘ 
eG (Tra) Sy we 
Q n3/2 
=> 


CbJSo (T+ 2h/b)"? 


Q 2/3 
h = | —— 1 + 2h/b)1/3 
> (ais) (1 + 2h/b) 


Here, 
h = 0.4028(1 + 0.9091h)?/9 

Iterate (from, e.g., h = 0.4028) to get normal depth 
h, = 0.4518 m 


For the critical depth, 


Coe = 2.045 m? s71 
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Then, 


2\ 1/3 
c= (5) = 0.7526 m 
g 


The normal depth is supercritical since hp < he (i.e. shallower and faster than the case Fr = 1, 
and hence its Froude number must be > 1). Hence, by definition, the slope is “steep” at this 
discharge. 


Answer: normal depth = 0.452 m; critical depth = 0.753 m. 


(b) First find the depth just upstream of the sluice gate by requiring the same head on both 
sides: 


Vy Vz 
Zsı + 2g = Z52 + 2g 
Q? Q? 
h +— >> =h, + 
1 2gb?h2? ° 2gb?h2 


With hz = 0.35 m: 
0.2132 
hy + —~=— = 2.090 
hy 
Rearrange for iteration to the subcritical (deep) solution: 
0.2132 
Iterating from, e.g., h4 = 2.090 gives h, = 2.039 m. 


The depth far upstream is supercritical, that upstream of the sluice is subcritical. Hence, there 
must be a hydraulic jump upstream, followed by GVE up to the gate. 


normal 


normal 


Answer: 2.04 m. 


(c) Upstream, since the channel is long, the jump goes from normal depth to the corresponding 
sequent depth in the hydraulic jump: 


hy = 0.4518 m 
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= 4,527 m s71 


Fr, = = 2.150 


T 
neat 


The GVF is subcritical, so work upstream from họ = 2.039 m just upstream of the sluice to 
hz = 1.166 m just downstream of the jump, using 2 steps in the GVF equation. (Note the 
renumbering of depths for this part.) 


TF JIF8FrA) = 1.166 m 


1.166 — 2.039 
= ——————_ = —0.4365 m 
2 
With lengths in m throughout: 
v2 Q? 0.4263 
Fr2 = — = ——————— = ——_ 
E gb?h? h? 
2 (1+ 2h/b) E _4 (4 + 0.9091h) 
Hence, 
0.4263 
dx 1=Fr 1- 
dh So- 5+ 0.01 — 6.537 x 10-4 EE 
and 
A (=) Ah 
x = — 
dh mid 
Working is set out in the following table. 
i hi Xi hmia (dx/dh)mia Ax 
2.039 0 
1.821 95.69 —41.77 
1 1.6025 —41.77 
1.384 88.87 —38.79 
2 1.166 -80.56 


Answer: distance to the (upstream) hydraulic jump = 80.6 m. 
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Section 4.2 


Example. (From White, 2006) 

A pencil point piercing the surface of a wide rectangular channel flow creates a wedgelike 25° 
half-angle wave, as in the figure right. If the channel has a Manning’s n of 0.014 m7 s and 
the depth is 350 mm, determine: (a) the Froude number; (b) the critical depth; and (c) the 


critical slope. 


a = 25° 
n=0.014m71/3 s 
h=0.35m 


(a) 


Answer: 2.37. 


(b) 


Fr = —— 


Vgh 
=> V=Fr/gh = 2.366 x v9.81 x 0.35 = 4.384 m s™t 


= q=Vh =4384x035 =1.534m? s7! 
gx 1/3 

> h= (£) = 0.6213 m 
g 

Answer: 0.621 m. 


(c) The critical slope is that at which the normal depth is equal to the critical depth. 
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Normal flow: 


1 . 
q=Vh where Y=— Re! 3ç1/2 (Manning) R, „=h (“wide” channel) 
n 


> q= * h2/351/2h 
n 
Rearranging, and setting a depth h equal to the critical depth (he = 0.6213 m), 


_ (nq)? (0.014 x 1.534)? 


> S 


Answer: 2.25x10°. 
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Sediment Transport; Section 2.1 
Note: removed from the present curriculum because of the reduction in teaching hours. 


Example. 

An undershot sluice is placed in a channel with a horizontal bed covered by gravel with a 
median diameter of 5 cm and density 2650 kg m°. The flow rate is 4 m° s™! per metre width 
and initially the depth below the sluice is 0.5 m. Assuming a critical Shields parameter Tčņit of 
0.06 and friction coefficient cp of 0.01: 

(a) find the depth just upstream of the sluice and show that the bed there is stationary; 


(b) show that the bed below the sluice will erode and determine the depth of scour. 


d=0.05m 

ps = 2650 kgm™? (s = 2.65) 
q = 4m? s™t 

T% = 0.06 

cp = 0.01 


(a) The bed is mobile if and only if the bed shear stress exceeds the critical shear stress for 
incipient motion. (Alternatively, one could find the velocity at which incipent motion occurs.) 


First find the critical shear stress for incipient motion. The local shear stress can be found by 
finding velocities using open-channel flow theory and using the friction coefficient. 


The critical shear stress for incipient motion is 


Terit = Tet(bs —p)gd = 0.06 x (2650 — 1000) x 9.81 x 0.05 = 48.56 Pa 


In the accelerated flow just below the sluice gate, 


h = 0.5 m 
E ok 42 -1 
V = A O05 a 8ms 
The total head relative to the undisturbed bed is thus (since z, = h initially): 
V2 2 
H=z,+— =05+22 =3.762 
4s 39 2x 9.81 = 


Upstream of the sluice one must have the same total head: 
2 2 


q 
H = — — h — 
ZF 2g + 2gh? 
but now we look for the subcritical (deep) solution: 
2 
q 
h=H- 
2gh? 
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Substituting numerical values, in metre-second units: 


0.8155 
Iterating (from, e.g., h = 3.762) gives 
h = 3.703 m 


The velocity upstream of the sluice is 


> a 
Lm ~ 3.703 


and the bed shear stress upstream of the sluice is 


= 1.080 m s™t 


1 1 
Tp = Cf (5 ev?) = 0.01 x 3 x 1000 x 1.080? = 5.832 N m~? 
Since Tp < Terit the bed here is stationary. 


Answer: depth upstream of sluice = 3.70 m; as demonstrated, the bed here is stationary. 


(b) Initially, beneath the sluice, 


V=8 m s™t 
1 1 2 -2 
ty = q (32v?) = 0.01 x 5 x 1000 x 8 = 320 Nm 
Since Tp > Terit the bed here is initially mobile. 


The bed will continue to erode, the flow depth h increasing and velocity V decreasing until 


Tp = Terit 
At this point, 


1 2 
GPV ) = Terit 


whence 
2 Terit 2 48.56 
SA Sf ee eet 
C p 0.01 ` 1000 S 
ret ee ipg 
=y gae AE 


This is the depth of flow. Since the downstream water level is set (at 0.5 m above the 
undisturbed bed) by the gate the depth of the scour hole is 


1.284 — 0.5 = 0.784 m 


Answer: 0.784 m. 
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